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Loop condensed phases are scale-invariant quantum liquid phases ol matter. These phases include 
topologically ordered liquid phases such as the toric code as well as critical liquids such as the 
Rokhsar-Kivelson point of the quantum dimer model on the square lattice. To investigate the exent 
to which non-local geometric observables capture a signature of the non-local quantum order present 
in these phases, we compute geometric properties of such loop condensed states using directed loop 
Monte Carlo calculations. In particular, we investigate the loop condensed nature of ground state 
of the square lattice quantum dimer model at the Rokhsar-Kivelson point and compare with other 
loop condensed states on the square lattice, including those of the toric code and fully packed 
loop model. The common features of such liquids are a scale invariant distribution of loops and a 
fractal dimensionality of spanning loops. We find that the fractal dimension of the loop condensate 
of the square lattice quantum dimer model at the Rokhsar-Kivelson point is 3/2, which provides 
quantitative confirmation of the effective height model that is commonly used to describe this critical 
dimer liquid. 
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I. INTRODUCTION 

The discovery of quantum liquid phases of matter 
that possess quantum order without breaking conven- 
tional symmetries has demonstrated that ordered quan- 
tum phases exist outside of Landau's paradigm of con- 
ventional phases of matter. To understand the formation 
of order in such quantum liquid phases, the mechanism of 
and loop and strinMiet condensation has been proposed 
by several authors 1"^. In this description of quantum 
liquids, the effective degrees of freedom are extended ob- 
jects (loops or string-nets) that fluctuate on all length 
scales and generate the non-local quantum entanglement 
present in such phases. In particular, loop condensates 
can describe Z2 topologically ordered phases as well as 
critical quantum liquids^. 

Since these quantum liquids cannot be described by 
a local order parameter, the presence of order in these 
phases must be characterized by other quantities. For 
example, signatures of topological order appear in the bi- 
partite entanglement entropy of a topologically ordered 
quantum liquid^ ^. A relevant question is whether these 
quantum orders have signatures in more conventional ob- 
servables. To make direct contact with the loop conden- 
sate picture of a quantum liquid, we seek to understand 
to what extent quantum liquid phases can be charac- 
terized directly in terms of the effective loop degrees of 
freedom. For example, in addition to local loop correla- 
tion functions, loop condensates can be described by the 
fractal dimension of the loop ga^. 

Lattice models that generate local constraints on the 
ground state often possess loop condensed liquid phases. 
The toric codc^ ground state can be described as having 
a local constraint requiring tha t an even number of loop 
segments touch each vertejP'^. The Z2 topologically or- 
dered ground state is simply the loop condensed phase 



of this (intersecting) closed loop subspace. Other models 
with local constraints may not explicitly be loop models, 
but the local constraint often may be readily mapped to 
a clos ed loo p constraint. For example, quantum dimer 
model^i^*^ impose a hard constraint at each vertex re- 
quiring a single dimer to touch each vertex: however by 
superimposing dimer configurations on a reference dimer 
configuration, a dimerization can be mapped onto a loop 
covering of the lattic^^^^^. On the square lattice, the 
quantum dimer model has an isolated critical dimer liq- 
uid poinlJSnnill^ but on the triangular lattice there is a 
gapped topologically ordered liquid phas^^SHU, Jn this 
work and in Ref. |23] we show that these dimer liquid 
phases can be viewed as loop condensed phase^^H, 

In this paper we investigate loop condensed states that 
arise in locally constrained models on the square lattice. 
We explicitly demonstrate the loop condensed nature of 
several liquid states in geometrically constrained models 
on the square lattice, including the square lattice quan- 
tum dimer model. In particular, we compute the frac- 
tal dimension of the underlying loop gas, and examine 
the fractal dimension as distinguishing feature of these 
phases. We find that the fractal dimension of the critical 
dimer liquid at the Rokhsar-Kivelson (RK) point of the 
square lattice dimer model is 3/2. This agrees with that 
of the contour loops of the effective height model descrip- 
tion of the RK point^^, providing quantitative validation 
of the height model. 



II. INTRODUCTION TO LOOP CONDENSATES 

A loop condensat^SEl jg g, scale invariant liquid state 
with fluctuating loops on all length scales. We may de- 
scribe such a loop liquid with both conventional local cor- 
relation functions and non-local geometric observables. 
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Local correlation functions include the loop-loop corre- 
lation function, C(r), and the one- loop correlation func- 
tion, Gi(r). The loop-loop correlation function C(r) de- 
termines the probability that loop segments exist at two 
points separated by r. This distinguishes between purely 
disordered loop liquids with exponentially decaying lo- 
cal correlations, and those with quasi-long range order 
due to local correlations that decay as a power law. The 
one-loop correlation function determines the correlations 
between defects that live at the end of a broken loop 
separated by a distance r. While these local correlation 
functions display the disordered nature of loop liquids, 
they do not capture the nature of non-local quantum or- 
der present in these phases. 

Non-local, geometric properties of loop condensates are 
described by the loop distribution P(s), and the fractal 
dimension D f , which are related to the two-loop correla- 
tion function 6*2 (r). Scale-invariance demands that the 
loop distribution, P{s), the probability that a loop has 
length s, scales as a power law of length: 

P{s)^s-\ (1) 

Additionally, the scaling of the length of a loop £ with its 
radius R{£) determines the fractal dimension Df oi the 
loop condensate: 

s{t)^R{£f^. (2) 

These geometric exponents may be related by a scaling 
relation to the two- loop correlation function G2{r) which 
determines the probability that two points separated by 
a distance r lie on the same loop'^^. Since the two-loop 
correlation function is fundamentally non-local (since it is 
determined by the existence of a single loop connecting 
the points), G2{r) may decay as a power law even in 
the absence of local correlations. The exponent X2 that 
governs the power-law decay G2 {r) ~ r~^^^ determines 
both the geometric exponents r and Dp^^. 

Df = 2-X2, T = l + —^. (3) 

2- X2 

In the remainder of this paper we will investigate to what 
extent non-local, geometric loop observables such as the 
loop distribution and the fractal dimension display sig- 
natures of the quantum order of loop condensates on the 
square lattice. 

III. LOOP CONDENSATION FROM LOCAL 
CONSTRAINTS 

Lattice models with local constraints can often be 
mapped to loop models, and consequently display a va- 
riety of loop condensed phases. Consider a lattice model 
with Ising degrees of freedom that live of the links of the 
lattice: the two states of each Ising degree of freedom may 
be considered to correspond to a link that is either oc- 
cupied or unoccupied by a dimer. A link / (un)occupied 
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FIG. 1. Examples of configurations with the local constraints 
on the square lattice. Clockwise from upper left: {Ci}, {C2}, 
{Co}, {Ce}, which have local constraints that restrict config- 
urations to have a one, two, an odd and an even number of 
dimers touching each vertex, respectively. 



by a dimer is the +1(0) eigenstate of the dimer num- 
ber operator ni = d^d^, with d^^ the dimer creation 
and annihilation operators. The simplest local constraint 
is to fix the dimer number at each vertex to be a con- 
stant no; we define the configuration space {Cn„} to be 
the set of configurations with no dimers touching each 
vertex. On a square lattice there are only two distinct 
non-trivial cases: rig = 1 (or equivalently rtp — 3) and 
uq — 2. {Ci} comprises fully-packed, hard-core dimer- 
izations of the lattice. {C2} comprises fully-packed, non- 
intersecting loop configurations with two dimers touching 
each vertex. 

All fully-packed, hard-core dimerizations in {Ci} can 
be mapped to a non-intersecting closed loop configura- 
tion by choosing a reference dimerization Rq (see figure 
When a dimerization is superimposed over Rq , ev- 
ery vertex will be touched by one physical dimer and by 
one reference dimer. With the exception of links where 
a physical dimer coincides with a reference dimer, the 
combined physical and reference dimers form a closed 
loop covering of the lattice, where the loops comprise 
sequences of links that are alternately occupied by phys- 
ical and reference dimers, respectively. We may choose 
to define links with overlapping physical and reference 
dimers to be unoccupied in the loop picture; with this 
choice the physical dimer configuration that completely 
coincides with Rq corresponds to the empty loop config- 
uration. 

The simplest liquid state in these constrained Hilbert 
spaces is given by a generalization of the Rokhsar- 
Kivelson wavefunction of the quantum dimer modeP^. 
This generalized Rokhsar-Kivelson wavefunction I^'J^k) 
is an equal super-position of all configurations obeying 
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the constraint: 



I*r"k) 



|C„o) 



(4) 



On a surface with trivial topology, the sum in Q is 
over all configurations C„q that obey the given local con- 
straint and is the number of such configurations. For 
no = 1, I^rk) i^ Rokhsar-Kivelson state of the quan- 
tum dimer model. Diagonal expectation values of I^S*^) 
are equal to those of the corresponding classical modeP^. 
Consequently, the behavior of correlation functions and 
local order parameters of |^rk) determined by the 
statistical mechanics of the related classical models. Di- 
agonal correlation functions of |^rk) determined by 
the classical dimer model, and those of |^rk) by the clas- 
sical fully-packed loop model, with a loop fugacity equal 
to 1 (which maps to the 6-vertex model) ^^'^^ . 

Rokhsar and Kivelson defined a local Hamiltonian for 
which I^rk) *be exact zero-energy ground state at one 
point in the parameter spac^^^ 

i?RK = -t ^ (|D) (□! + h.c) +vY,{p) (□! + p) (□ 
p p 

(5) 

In Eq. ([5]) the sum is over all plaquettes p of the lattice. 
The "f term fiips dimers around a "flippable" plaquette 
(i.e. a plaquette with alternating occupied and unoccu- 
pied links) and the "u" term is a potential energy for such 
fiippable plaquettes. For v — t, the so-called RK point, 
I^rk) energy ground state of H^k- We note 

that 77rk is specifically constructed such that it con- 
serves the dimer number at a given vertex, and therefore 
does not violate the local dimer number constraint. 

On a surface of nontrivial topology, such as a torus, 
the local dynamics of -//rk will break {C„o} into distinct 
topological sectors. i?RK will not connect configurations 
in different topological sectors and thus generates a topo- 
logical degeneracy at the RK point. Indeed, a distinct 
ground state I^Prk) ^^"^ ^e defined within each topologi- 
cal sector by limiting the sum in Eq. Q to configurations 
within the topological sector. The nature of this topologi- 
cal degeneracy is manifest in the loop condensate picture: 
each topological sector corresponds to a distinct winding 
sector defined by the number of loops winding around 
each topologically nontrivial cycle of the surfac^^J, 

Moessner et al.^ have shown that a local dimer num- 
ber constraint generates a local U{1) invariance. If we 
define the operator n^ = where n/ is the dimer 

number operator on the link I, then all states \ipno) that 
obey the local constraint are invariant under the follow- 
ing local gauge transformation: 
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J}» (a) = exp (^ia (n„ - no)) , G^J" (a) \^no) = \^no) 



(6) 



Consequently, i?RK bas bee n sho wn to map to a U{1) 
gauge theory for both no = llsHUni and no = 



We now consider relaxing the dimer number constraint 
to a local dimer parity constraint. Here {Ge} and {Go}, 
which comprise configurations with a fixed dimer num- 
ber parity at each vertex, have an even or odd number 
of dimers touching each vertex, respectively. The phys- 
ical states \tpe,o) that are formed from superpositions of 
{IC'e.o}} Siie invariant under the gauge transformations 

GS =exp(±z™,), G^|V'e) = |^e) (7) 
G° = exp (iiTT (n„ -fl)) , G° IV^o) = IV-o) . (8) 

Note that the U{1) symmetry of G"" (a) has thereby been 
reduced to a Z2 symmetry in G^'"^^^. 

We may now define a Hamiltonian that commutes with 

Ge,o. 
„ ■ 



P l<£p 
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(9) 



where df = l/\/2(d;^ + d^) and d^' . This is the mag- 
netic term of the toric code HamiltoniarfS. In the toric 
code Hamiltonian the dimer parity constraint is imposed 
by a local energy cost at each vertex; here we take that 
energy cost to be infinite such that this is a hard con- 
straint. The ground state of Eq. ^ is the equal super- 
position of all configurations in the given parity sector: 
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TC/ 



(10) 



In the language of the Moessner et aL'^^-', l^'xc) ^^e 
even (odd) parity sector corresponds to the even (odd) 
Ising gauge theory; l^'f-^) is simply the toric code ground 
state. 

We may now understand the relationship between the 
U{1) gauge theory describing Eq. Q with a local dimer 
number constraint, and the Z2 gauge theory that de- 
scribes Eq. (10 1 with a fixed dimer parity constraint. 



Since the square lattice is bipartite, the links can be ori- 
ented to always point from one sublattice to the other; 
dimers can be viewed as carrying a flux into or out of 
a vertex according to this orientation. The riy = 2 con- 
straint acts as Gauss's law, such that there is a flux of +2 
and —2 on each vertex of sublattice A and B respectively; 
correspondingly we can consider there to be a static back- 
ground charge of ±2 on the two sublattices. Now consider 
if I'I'rk) is doped with dynamic n„ = 0,4 vertices; these 
will act as charge ±2 objects, depending on the sublat- 
tice. Allowing n„ = 0,4 vertices will transform I^'rk) 
into which is the toric code ground state, and is 

known to be described by a Z2 gauge theorj^ . 

This transition from a [/(I) gauge theory to a Z2 gauge 
theory via the introduction of charge 2 objects follows 
the well know prescription of Fradkin and Shenker, who 
showed that coupling a U{1) gauge field to charge > 1 
matter field can reduce the U{1) gauge symmetry down 
to Zf^^. The same picture applies in the odd parity 
sector, where the QDM can be viewed as a U{1) gauge 
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theory with ±1 background charges on the sublattices. 
In this case, introducing n^, = 3 vertices is equivalent 
to allowing charge 2 matter fields and generates l^'xc)' 
which is described by the odd Ising gauge theory. This 
construction leads one naturally to ask whether other 
Zjv models may live at the RK points of geometrically 
constrained models on the square lattice. In particular, 
doping the QDM with charge ±3, riy = A vertices might 
lead to a liquid phase, purely as a result of geometrical 
constraints. 



IV. SQUARE LATTICE QUANTUM DIMER 
MODEL 

The ground state of the RK point of the square lat- 
tice QDM is a gapless critical dimer liquid, with power- 
law decaying dimer-dimer correlationgii^. As described 
above, a dimer model may be mapped to a closed loop 
model by introducing a reference dimerization. As such, 
we may view the RK point of the square lattice QDM 
as a loop condensed liquid phase. There is a well known 
phenomenological height model description of the square 
lattice RK wave functiorP^H2l!. Here we will give a vari- 
ant of this height model picture that makes the meaning 
of the loop picture more transparent (see figure [2| : 

• We first we orient the links such that the arrows 
point from sublattice A to B and choose a reference 
plaquette to assign a height of /i = 0, as in figure 

13 

• Heights h{r) are assigned to all other plaquettes by 
starting with the reference plaquette and moving 
along any path crossing the links, where we assign 
a Sh for each link, with \Sh\ = 1. 

• If a link crossed with the arrow pointing to the right 
is occupied by a physical dimer, Sh = +1, occupied 
by a reference dimer 5h = —1, and vice versa if the 
link points to the left. 

• Flipping parallel dimers on a plaquette will change 
the local height by ±1. 

With this mapping, the transition loops formed from the 
combination of the reference and physical dimerizations 
act as the contour loops of the height field, since the 
height only changes when a transition loop is crossed (see 
Fig. [2]). 

We then propose an effective gaussian action S'rk for a 
coarse grained height field which captures the local height 
fluctuations and describes the dimer liquid at the RK 
point: 



Srk = / dr 



\^h\' + V{h). 



(11) 




The statistical weight for a given configuration of heights 
is then proportional to ea;p(— S'rk)- In Eq. ( 11 ), K is the 



FIG. 2. An example of the mapping of the square lattice 
QDM to a height model. The green dimers represent the phys- 
ical dimers and the blue dimers are the background dimeriza- 
tion. The lower left plaquette is chosen to have ho = Q and 
all other heights are determined relative to this, by following 
a path from ho through the links of the lattice. If a dimer is 
crossed with the arrow pointing to the right (left) the height 
changes by +1 (—1) and vice versa for a reference dimer. Flip- 
ping dimers on a flippable plaquette changes the local eight 
of the plaqeutte by ±1. The transition loops formed by al- 
ternating green and blue dimers are the contour loops of the 
height field. 



stiffness and the first term captures the fiuctuations of 
the height field. V{h) is a locking potential which favors 
certain ordered height configurations. At the RK point 
on the square lattice, the locking potential is irrelevant 
in the renormalization group sense, and therefore the ef- 
fective action is gaussian. This gaussian action describes 
a rough phase of the height model, where h fluctuates 
along the lattice. Such a phase can be shown to have 
power law correlation functiond^MH! Consequently, the 
height model can capture the critical correlations of the 
dimer liquid at the RK point, with an appropriate choice 
of the stiffness K. We note that this height model de- 
scription is purely phenomenological-currently there is 
no well known microscopic derivation of this action for 
the QDM. 

The transition loops to the reference dimerization take 
on a special meaning in the height model. In particular 
they are the contour loops of the height model, in that 
the height only changes when a transition loop is crossed 
(see Fig. [2|. The fractal dimension of the contour loops 
of a gaussian height model has been shown to be univer- 
sal (independent of the stiffness K) and equal to 3/2 by 
Kondev and Henley^*^. 

We have computed the distribution of the transition 
loops of the square lattice QDM at the RK point with 
directed loop Monte Carlo calculations'^^ . We first con- 
sider the computed loop distribution shown in Fig. [3j We 
see that P{s) displays a power law behavior over loops 
on length scales s « L^, which is indicative of the loop 
condensed nature of the dimer liquid at the RK point. 
The best fit power law gives r = 2.32 ± 0.02. In the 
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FIG. 3. The loop distribution function P (s) for the square 
lattice dimer model at the RK point for a lattice of linear 
dimension L = 512. P{s) displays a clear power law over 
length scales s <^ . The line shows the best fit power law 
with T = 2.32 ± 0.02. 




FIG. 4. The distribution of the longest loops , P{sm) for 
the square lattice quantum dimer model at the RK point on 
an L = 64 lattice. The broad distribution over many length 
scales is indicative of the loop condensed nature of the state. 
The location of the peak scales with system size as shown in 
figure |5] 



effective height model, X2 = 1/2 and (from the scaling 
relation given in Eq. ^) r ~ 7/3, so our result is quan- 
titatively consistent with the predictions of height model 
description. 

For loops approaching a length L^, the longest length 
scale of the lattice, we see a clear deviation from the 
power law behavior. Fig. [4] shows the distribution of the 
longest loop; this displays the broad distribution over 
many length scales that is characteristic of a loop con- 
densate. The location of the peak of P{sm) scales with 
system size, but not with an integral power (e.g. 
or L^) which would be expected for dilute loop or sym- 
metry broken phases, and therefore these spanning loops 
have a fractal dimensionality. Fig. [5] shows the finite size 
scaling of the fractal loop: the best fit power law gives 
Df = 1.502 -I- ±0.002 which is clearly consistent with 
the universal height model fractal dimension prediction 
of 3/2. 

These calculated geometric exponents provide quanti- 
tative confirmations of the height model description of 
the square lattice QDM. In contrast with the local corre- 
lation functions (for which the value of K must be chosen 
to match the form of the height model correlations to the 
computed dimer model correlations functions) the uni- 
versality of the fractal dimension of contour loops means 
that this quantitative agreement does not require tuning 
any parameters. Additionally we note that the fractal 
dimension of the transition loop distinguishes the critical 
liquid RK point from the gapped liquid of the triangular 
QDM, for which we have determined the fractal dimen- 
sion to be 7/^. 
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FIG. 5. Finite size scaling of the length of longest loop, sm 
of the ground state of the RK point of square lattice QDM 
(|*?rc))- The best power law fit gives Df = 1.502 ± 0.002 
which is consistent with the fractal dimension of contour loops 
of a gaussian height model. 



V. OTHER LOOP CONDENSATES ON THE 
SQUARE LATTICE 

A. Toric code 

The ground state of the toric code, l^'xc)' '^^Y be de- 
scribed as a loop condensate on the square lattic^. This 
gapped liquid phase has exponentially decaying spin-spin 
correlation functions. The ground state subspace is that 
of closed loop coverings of the square lattice, with the ex- 
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FIG. 6. Finite size scaling of the longest loop in toric code 
ground state (|*tc))- We find D/ = 1.7502 ± 0.0002, con- 
sistent with the universality of the fractal dimension of the 
gapped Z2 liquid. 



ception that two loops may meet at a vertex. Since l^'xc) 
is an equal superposition of all loop coverings, we may 
relate it to the 0(1) loop model, by choosing a resolution 
of the four-loop vertices. Here we choose an orientation 
of the loops, such that loops do not cross at an n„ = 4 
vertex (see figure [I]). In figure |6j the finite size scaling of 
the largest loop of the toric code wave function is plot- 
ted. From this we have extracted the fractal dimension to 
be Df = 1.7502 ± 0.0002, which agrees with the known 
value for the 0(1) model, Df = 7/42S1. We have also 
computed the fractal dimension for the odd Ising gauge 
theory described by and found that it is also 7/4. 

Previous work has shown that the fractal dimension of 
Z2 topological phases of the triangular lattice quantum 
dimer model and the honeycomb loop models both have 
fractal dimensions of 7/4^ This suggests that the frac- 
tal dimension of 7/4 is universal for the Z2 topological 
liquid phases. 



B. Fully packed loop model 

The fully packed loop model is a critical liquid state 
with algebraically decaying loop-loop correlation func- 
tions^^ Just as with the RK point of the square lat- 
tice QDM there is also a phenomenological mapping of 
this model, which describes the RK wave function with 
a 2-dimer constraint at each vertex, I^'rk); a height 
modei^'' . Similar to the dimer model case, the loops 
are oriented and act as the contour loops of the height 
field. However, in this case the effective action S'fplm is 
a Gaussian action augmented by an additional term S' 
that is relevant (in the renormalization group sense), and 
which acts as a background charge in the Coulomb gas 
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FIG. 7. Finite size scaling of the length of the longest loop of 
the fully-packed loop model that describes the 2-dimer con- 
strained RK wave function, |^rk)- The best power law fit 
gives Df — 1.7501 ± 0.0002 which is consistent with the frac- 
tal dimension of the effective height model. 

description of the height modeP^^^^l^ 

^FPLM = J dr^(^^\^hf + V{h)') +S'. 

In this situation the universality of the Df = 3/2 for 
a pure Gaussian model does not apply. Indeed, for the 
fully packed loop model, the height model mapping pre- 
dicts a fractal dimension of = 7/4. We have con- 
firmed this with directed loop Monte Carlo calculations 
of I^'rk) summarized in Fig. [Tj a power law fit yields 
Df = 1.7501 ± 0.0002. While the fractal dimension does 
not distinguish the fully packed loop condensate from 
the gapped Z2 loop condensed phase, these two phases 
are easily distinguished by the density-density correlation 
function, C{r). In particular, the density-density corre- 
lation function displays a power law in the critical fully 
packed loop model, while it vanishes exponentially in the 
gapped liquid phase. 

VI. CONCLUSIONS 

We have studied the geometric properties of the tran- 
sition loops of ground state of the RK point of the square 
lattice QDM. The power law distribution of loop lengths 
as well as the fractal dimensionality of loops that span the 
finite system indicate that this critical dimer liquid state 
is fundamentally described by a loop condensate. Our 
numerically computed geometric exponents agree with 
the predictions of the standard height model description, 
which gives a quantitative confirmation of the relevance 
of this phenomenological action. We have compared the 
loop condensate of the square lattice QDM with those of 
the toric code and fully packed loop models on the square 
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lattice. The fractal dimension of 3/2 of the square lattice 
QDM distinguishes it from the gapped Zi topologically 
ordered phase of the toric code on the square lattice and 
that of the QDM on the triangular lattice, both of which 
have fractal dimension of 7/4. However the critical loop 
condensate of the fully-packed loop model has the same 
fractal dimension as these gapped phases, so the cur- 
rent work demonstrates that these geometric exponents 



may not always distinguish loop condensed phases. Sim- 
ilar analyses of other quantum liquids phases may prove 
fruitful for characterizing other loop condensed phases. 
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